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We propose a new semiclassical approach based on the dynamical mean field theory to treat the 
interactions of electrons with local lattice fluctuations. In this approach the classical (static) phonon 
modes are treated exactly whereas the quantum (dynamical) modes are expanded to second order 
and give rise to an effective semiclassical potential. We determine the limits of validity of the 
approximation, and demonstrate its usefulness by calculating the temperature dependent resistivity 
in the Fermi liquid to polaron crossover regime (leading to 'saturation behavior') and also isotope 
effects on electronic properties including the spectral function, resistivity, and optical conductivity, 
problems beyond the scope of conventional diagrammatic perturbation theories. 
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I. INTRODUCTION 

The present-day understanding of electron-lattice in- 
teractions in metals is based on the assumption that typ- 
ical phonon frequencies are small relative to typical 
electron energies t, so that electronic properties may be 
calculated as an expansion in the adiabatic parameter 
7 — Q./t. Tiiis was first exploited by Migdal and Eliash- 
berg (MEp to derive a self-consistent set of equations 
for electron and phonon self-energies. Since ME theory 
neglects order 7 vertex corrections it wrongly predicts a 
zero result for most isotope effects on electronic proft. 
erties (notable exceptions being the superconductinga 
and charge-orderingj transition temperatures). Further- 
more, ME theory assumes that the underlying electronic 
groundstate can be described by Fermi liquid theory. 
However, if the electron-lattice coupling strength X = A/t 
exceeds a critical value (of order 1) the conduction elec- 
trons are believed to form "small polarons" so that the 
electronic groundstate is fundamentally reconstructed. 
In this case, although an expansion in 7 exists, the start- 
ing point is not clear. 

Signatures of polaronic behavior have been ob- 
served in certain complex narrow band materials. 
Well-known examples are phonon mediated supercon- 
ductors with high traxisition temperatures, likep.BiO 
based superconductoijS,Q alkali-doped Cgo fullerides,l3 and 
the A15 compounck,Ll and 'colossal magnetoresistance' 
(CMR) manganitesjj Despite much work many questions 
remain open: as examples we cite a possible 'saturation' 
of the resistivity at high temperatures ,0 unexpected mid- 



infrared transitions in the optical conductivity,!^! and 
large isotope effects on electronic properties. For in- 
stance, the CMR manganite (Lao. 25Pro. 75)0. yCao.sMnOa 
undergoes a metal-insulator transition upon ^^O - ^^O 
isotope substitution.EI A calculational technique which 
can address with these issues is urgently needed. 

The interaction of conduction electrons with local lat- 
tice fluctuations may be d£scribed by the dynamical 
mean field theory (DMFT).E£l In this approximation the 
electron self energy is momentum independent. Since 
the full frequency dependence is retained the correspond- 
ing one-electron Green function shows both coherent and 
incoherent features. Recent applications of DMFT to 
the electron-phonon problem include a systematic expan- 
sion in nawjei:a of 7,0 several weak-coupling expansion 
schemes,0OEj and a "classical" (7 = 0) study of the 
Fermi liquid to polaron crossover. (Treating the 7 = 
limit- without assuming DMFT has not been possible 
yet.li3) The expansion of Ref. ^ is restricted to very 
low temperatures and cannot access the intermediate or 
strong coupling regimes where important physics such 
as resistivity saturation or polaronic effects is expected. 
Because of the mismatch between phonon and electron 
energy scales numerical simulationalB are difficult to im- 
plement accurately, except in the "antiadiabatic" limit 
7 1 of unclear physical relevance. 

Here we propose a new "semiclassical" method for 
treating intermediate and strong couplings in the physi- 
cally relevant limit 7^1. We treat the classical (static) 
phonon modes exactly (as in Ref. |4|) and expand the 
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quantum (dynamical) modes to second order. The re- 
sulting action consists of a classical part plus an effective 
semiclassical potential. The method is valid for a wide 
range of electron-phonon couplings and resulting ground- 
states. Here we consider only electron-phonon interac- 
tions (i.e. the Holstein model). However, the method 
can be generalized (along the lines of Ref. [Tl| , [l6| ) to study 
phonon effects in the presence of strong electron-electron 
interactions. We focus on the crossover regime between 
Fermi liquid and polaron physics which was not easily 
tractable by previous methods. The remainder of the 
paper is organized as follows. In Sec. H we introduce the 
Holstein model and the semiclassical method. In Sec. 

III we discuss additional simplifications which are con- 
venient for numerically implementing the method. Sec. 

IV studies analytical limits and defines the observables 
(spectral function, resistivity and optical conductivity) 
we consider. Numerical results are presented and dis- 
cussed in Sec. V. We conclude in Sec. VI. 



II. MODEL AND FORMALISM 

In this paper we study the Holstein modelEl of elec- 
trons interacting with local lattice fluctuations: i/hoi = 
ffci + Hpi, + i?oi-ph with 



(2) 

Hel-ph {4^C^^ - Tlj , (3) 



where we have absorbed the electron phonon coupling 
into the phonon coordinate r which thus has dimension 
of energy as does the phonon stiffness parameter A. In 
our numerical calculations we further specialize to a mean 
electron density per spin direction of n — 1/2 (this im- 
plies /U = 0) and we use the semicircular density of states 
p{e) = {l/N)j:j:S{e- ej^ = l/{27Tt^)VW^ (per 
spin), where is the Fourier transform ofU-j. However, 
our approach is not restricted to these cases. The possi- 
bility of charge-density wave formation (due to a nested 
Fermi surface) will not be considered in this work, but 
could also be considered via our method. 

The fundamental assumption of DMFT is the momen- 
tum independence of the electron self energy S(p, i LUn) — > 
E„ = I](zw„). This implies that the physics can be de- 
rived from an impurity model specified by the action 

S{{rk}; {cn}) = (T/2A) 5] rfe (1 + ul/Q^) rl~nr^ 

k 

-(2s + 1) Tr In [c„ 5nni - T r^-m] ■ (4) 



The Cn are mean field functions which describe the con- 
duction electrons and depend on odd Matsubara frequen- 
cies Ldn — 2ttT {n + 1/2). The rfc are bosonic fields which 
describe the phonons and depend on even Matsubara fre- 
quencies LUk = 2'KTk. The factor (2s -I- 1) is the spin 
degeneracy. The partition function may be written as a 
functional integral over the bosonic fields 



J V[r] exp(-S'). 



(5) 



It is a functional of the mean field function c. The im- 
purity Green function Q and self energy S are defined 
by 



Gn = 



1 SlnZ 



1 



2s +1 SCn C„-i;„({c„}) 



(6) 



The mean field function is fixed by equating the local 
Green function (Gj-*"-(^) = / (0 ZJ^T^^) of the 
original lattice model to Gn- For a semicircular density 
of states one findstj 



iuJn+ IJ.- Gn ({C„}) 



(7) 



The effective action (^ defines an interacting quantum 
field theory which cannot be solved exactly. Here we 
expand the action to quadratic order in the quantum 
modes LOk {k ^ 0) while treating the classical modes ojq 
exactly: 



Tr In [c„ Snm 

n^rn 



Trn-rn] ~ Tr In [c„ - T ro] 



Tro) 



' n—m 



m — n • 



(8) 



The effective action S ^ Sq + J2k>o separates into 
a classical part 5*0 and quadratic quantum parts Sk>o, 
both of which depend on rg: 



So 
Sk>o 
where 

and 



T 2 
2A 
T 



n ro - (2s + 1) Tr In [c„ - T vq] 



Tk [ujl + nl{ro)\ rl, 



nl[ro) = n' [i + nfe(ro)] 



Hfe(ro) = (2s + l)Ar^. 
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1 



Tro Cn-k - Tro 



(9) 
(10) 

(11) 
(12) 



Since the dynamic action (10) is quadratic the k ^ 
modes can be integrated out exactly. The resulting par- 
tition function 



drP(r) 



(13) 
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may be written as a one-dimensional integral over the 
rescaled classical coordinate r = T ro and 



P(r) = exp|-5o(r)-^ln 

I fc>0 



(14) 



is the (unnormalized) probability that the classical co- 
ordinate takes the value r. Note that P{r) depends on 
the mean field parameters Cn which must be computed 
self-consistently. Performing the derivative with respect 
to Cn yields the local Green function 



1 



drP(r) 



1 



1 



(15) 



The function Fn is given by 



Fnir)=ATj2 



wl+Cll{r) \Cn-k-r Cn+k 



(16) 



Eqs. ([7|) and (15) form a complete set of equations which, 
in principle, can be solved for the local Green function 
on the imaginary Matsubara axis. However, the direct 
numerical implementation is difficult, especially if the lo- 
cal Green function is required for real frequencies. In 
the following section we discuss physical approximations 
which considerably simplify the numerical task and allow 
for analytical insights. 



III. NUMERICAL IMPLEMENTATION 

A. Static approximation 

The rapid increase of Sk>o as the phonon frequency 
is increased above the Debye frequency SI, combined 
with the observation that life varies with frequency on the 
scale i, means that when computing the c„ [via Eq. (|^)] 
as well as static and energetic quantities we may replace 
nfe(ro) by its zero frequency limit no(ro). (For a more 
formal justification see Refs. 0,ll|.) Then Sk is the ac- 
tion of a harmonic oscillator with renormalized frequency 
n{r) = floir) = n[l + U{r,uj = 0)]^/^ and the phonon 
distribution function reads 



static 



(r) = 



n{r)/2T 



cxp (-S'o(r)) 



sinh(r2(r)/2T) 
The function F introduced above is given by 



Fn{r) 



A 02 



1 



2 n{r) 



coth 



2T 



2T 



nir) 



(17) 



(18) 



in the static approximation. Note that Fn ^ 1/T as 
T ^ oo. The form (p^ is used in the iteration of the 



DMFT equation (^ on the Matsubara axis. We use up 
to 512 positive Matsubara frequencies and treat high fre- 
quencies analytically. At the end of this first DMFT cycle 
we obtain the phonon distribution function P{r). 

The scheme is stable (the 'fc ^ 0' integration exists 
and P{r) > 0) provided n'^{r) > -iir^T^, in other 
words for all real renormalized phonon frequencies and a 
(temperature-dependent) range of imaginary ones. Note 
that ri2(r) increases as |r| is increased from zero. An 
imaginary phonon frequency corresponds to an instabil- 
ity of the ground state, but if T > |ri|/27r in this for- 
malism thermal fluctuations re-stabilize the ground state. 
For iniaginary ^(r) one must analytically continue Eqs. 

(HI), im to 



-Pstatic (?■) 



I»WI/2T 
sm{\n{r)\/2T) 



exp (-S'o(r)) 



(19) 



and 



Fn{r) 



A 



2 \n{r)\ Cn 



2T 



cot 



Mr)\ 
2T 



(20) 



respectively. The form (|l9|) has to be used for low T, 
large A, and small r (for example, the case A = 2.25, 
n = 0.1, and T < 0.2 [see below]). 

As will be discussed in more detail below an imag- 
inary renormalized phonon frequency is a precursor of 
polaronic effects. When the temperature is lowered we 
have |J7(r)|/2T tt. If we introduce a small deviation 
e(r) (increasing with r) via |f2(r)|/2T = tt — e(r) the 
semiclassical phonon distribution 



^'static(?') = -TT exp (-6*0) 
e(r) 



(21) 



is enhanced around r « in the polaronic phase. This 
may be interpreted as a precursor effect of a "pola- 
ronic band" where quantum oscillations allow for nearest 
neighbor polaron tunneling. An analysis of the regime 
T < |r2(r)|/27r requires a different treatment not given 
here. 



B. Dynamic corrections 

For dynamical quantities, like the optical conductivity, 
frequency corrections to IIq will be important at low tem- 
peratures (they lead e.g. to finite electron lifetimes) and 
the static approximation is not sufficient. Employing the 
full expression (|6|) in the calculation of the local Green 
function is cumbersome. To obtain Q on the real fre- 
quency axis a numerically challenging continuation from 
the imaginary frequency axis has to be performed. We 
thus prefer to do all our calculations on the real axis. We 
begin by expanding 
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dnlir) 



dujk 



(22) 



fc=0 



On the real frequency axis this corresponds to an expan- 
sion in small frequencies which dominate dynamic quan- 
tities at low temperatures. The second term in Eq. ( |2^ ) 
is imaginary and leads to a damping of the phonon modes 
and to small (relative order 7) corrections to physical re- 
sults. We will neglect this feedback of the electrons on 
the phonon system and keep only the frequency renor- 
malization, i.e. the first term in Eq. (p2|). This leads 
to 



and 

F„(r)=Ar^ 



Pir) = Pstatic(r) 



(23) 



fc>0 *- \ J \ n K 



Cn+k - r 



(24) 



Note the difference to the static approximation [Eq. (18)] 



where the finite frequency transfer of tuk in the effective 
electron fields c„_fc and Cn+k is missing. Eq. ( p^ ) defines 
the correct F for calculating dynamic quantities from 
Pstatic(?')- In the following we will restrict ourselves to 
the regime of stable quantum modes where Pstatic(^) is 
well defined down to T = 0. It is easy to analytically 
continue the function F to the real axis. The real and 
imaginary parts read 

, . A 

F (a;,r) = — — X 

h{Q.r) V'iuJ + Vlr)- bi-rir) V' [uj - Q.,.)- 



{2T/nr)V {uj) 

I 



-P I deV (e)/(e) 



and 



b{nr)r (to + Qr) -b{-nr)r (to-nr)- 

{2T/nr)V"iLu) - 



(26) 



Here jiuj) and biio) are the Fermi and Bose distributions, 
respectively. V and V denote the real and imaginary 
part of 7^(0; -I- «0+) = \l\c{uj) — r], respectively, fir = 
51 (r) and the integral in Eq. (ESI) is a principal value. If 
the propagators V exhibit particle-hole symmetry then 
F'(w) = -F {-w) and f"(w) = f" {-uj). The real part 
is an odd function of frequency and the modifications 
arising from a finite 17^ may be ignored at w « 0. To 
leading order in f2r we are left with 



(27) 



if{LU + nr)V {uj)-if{uJ-nr)V (uj) - {2T / ^r) V (uj)] 

Dynamic corrections reduce the imaginary part of F at 
small frequencies. For T — and r = the static ap- 
proximation yields a non-zero value F ~ 7 for small 
UJ whereas including finite frequency corrections give 
Pt=o = in a shell —Hq < uj < fto. 

The results of this section constitute an expansion in 
7, and have errors of order 7^. However, for temper- 
atures smaller then the renormalized phonon frequency 
the (quantum modified) thermal fluctuations become of 
order 7^ and one of the order 7^ errors induced by the 
semiclassical approach leads on the real axis to a self en- 
ergy whose imaginary part changes sign, so the self en- 
ergy becomes non-causal [see below Eq. (^)]. We thus 
suggest two complementary ad hoc schemes for iterating 
the DMFT equations on the real frequency axis. Both 
schemes result in a causal electron self-energy and 0(7^) 
or 0{'^T) corrections to physical quantities. 

(i) Setting F{uj, r) = for all uj and r. This approach 
gives the correct high temperature behavior and is also a 
good approximation for T = 0, and we expect it to give 
us an idea how the high and low temperature regimes 
of dynamic quantities are connected. Furthermore, the 
imaginary part of F does indeed vanish for small enough 
frequencies at zero temperature — as we would expect on 
physical grounds since the resistivity should drop to zero 
in a Fermi liquid state. 

(ii) The violation of causality is only of order 7^, as 
we will show in the next section. In order to obtain a 
well-defined electron self energy we may therefore add a 
heuristic impurity scattering rate r ~ 7 to the self energy: 



(25) 



Q^{g-^+iT) ' (28) 
We discuss a proper choice of r in the next section. 



Note that when crossing over from the Fermi liq- 
uid to the polaronic regime the renormalized frequency 
r2(r — 0) tends to zero. We believe that for A = Ac 
the ^F = 0' approach (which is then practically identical 
to the classical case) becomes exact for properties which 
involve only electrons at the Fermi energy. This would 
imply a square root like resistivity p ^ T^/^ down to 
T = at this singular point. Eil For coupling constants 
slightly below Ac the suggested approaches (i) and (ii) 
should give a reasonable picture down to lowest temper- 
atures. 

To summarize, the effect of quantum phonons to or- 
der 7 is to add a semiclassical potential to the classical 
action. Static and energetic quantities, including the dis- 
tribution function of onsite lattice distortions r, may be 
calculated using a simplified ('static') semiclassical po- 
tential depending only on a renormalized frequency fl{r). 
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The approach is vahd for all coupling strengths A but is 
restricted to temperatures T > \^l{r = 0)|/27r in the po- 
laronic regime A > Ac ~ 2.2 where 51 (r) becomes imagi- 
nary. Some additional numerical effort (compared to the 
classical casell^) is required due to the calculation of the 
renormalized frequencies. Care has to be taken when cal- 
culating dynamical quantities at low temperatures due to 
the presence of an additional quantum scattering term. 
We propose two schemes for incorporating dynamic cor- 
rections on top of the simplified semiclassical action. For 
coupling strengths A < Ac, i.e. in the Fermi liquid regime, 
the crossover between the high and low temperature be- 
havior of transport properties can be studied. 



IV. DISCUSSION AND OBSERVABLES 

A. Fermi liquid to polaron crossover 

When discussing the influence of local lattice fluctua- 
tions on electronic properties one has to distinguish be- 
tween two regimes. For small values of A the conduction 
electrons are weakly renormalized quasiparticles. This 
is the Fermi liquid regime. For large A the conduction 
electrons are self-trapped by their own local lattice dis- 
tortkms. This is the polaronic regime. In the classical 
caselij (7 = 0) the main features of the "polaronic transi- 
tion" may be illustrated by expanding to quadratic order 
around r = 0: 



So 



1 



2AT 



(i + Ar2) 



with 



Fc = i2s + l)Tj2 



(29) 



(30) 



For later use we define 



A, 



A 



classical 



1 + AF2' 



S^classical = ^ VT+AI^ . (31) 



As long as A < -l/r2 = A^^'''''''''''', 5*0 is minimized 
by r = 0. In the T —^ limit the dominant contribu- 
tions to Z come from small distortions. If A approaches 
^classical fj-Qj-^ bclow, the clastic energy vanishes and the 
system enters a polaronic state with (r) 0. The crit- 
ical Ac depends on band filling, temperature. A, as 
well as other parameters (such as electron-electron in- 
teraction strength). In the noninteracting case [where 
c° = ^ [uJn + \/uJn + 4*2] (upper half plane)] and aLT = 
we find Ac = 37ri/[(2s -h 1) 4]. It has been shownEJ that 
at high temperatures electrons are strongly affected by 
the presence of a local (polaronic) lattice instability. A 
strongly coupled system turns from a Fermi liquid into a 
polaronic insulator with activated conduction caused by 
the hopping of small polarons. In the crossover regime 
the system can be characterized as a "bad metal" . This 
intermediate state is the focus of the current paper. 



We now discuss the changes to this picture introduced 
by quantum effects. We first focus on the regime of sta- 
ble quantum modes, i.e. A < Ac. For small lattice dis- 
tortions we may expand the effective phonon frequency 
around r = 0. At half filling we find 



n^{r) (i-i-Ar2-f3Ar4r2 



(32) 



In the regime under consideration we can assume 1 -I- 
Ar2 > and use Eq. (17) to calculate the effective 
partition function. To second order in r and at low tem- 
peratures the semiclassical action reads 



S = 



1 



2AT 



A11V1 + AI^+ 1 + AF 



S{r = 0) + 



2 VTTAT^ 



2AT 



(33) 



The first (r = 0) term is the Migdal correction to the 
classical action given ioJilq. (^) . It arises from the one- 
phonon-loop diagramsO Quadratic and higher orders in 
the distortions have to be considered at higher temper- 
atures or stronger couplings, especially if the system is 
close to the instability Ac. In terms of the dimensionless 
parameters 



X^A/t, -f = n/t, 



(34) 



the leading corrections to the classical action are of the 
form A"i3^ggi^3^j 7ciassicai with n > 2, i.e. they are of order 
7ciassicai- Noticc howcvcr that the corrections can become 
arbitrary large if A ^ ^classic Jig 

The T — polaronic instability occurs when the coef- 
ficient of in Eq. ( ^3| ) vanishes, i.e. when (in rescaled 
units) 



1 + Ac {tT2 



Ag7 3(t3r4)/2 
^l-f Ac (tFa) 







(35) 



This equation has a solution at Ac = Ac (7) = —^2/1 + 
0{'^'^/^). Because r4 > (for the Holstein model at 
half filling r4 = (2s + 1) 8/(15 tt^^) in the noninter- 
acting limit) we see that quantum fluctuations increase 
the critical coupling needed for the polaronic instability 



by an amount 5\ 



classical 



2/3 



Note also that the sign of 



ri2(r = 0) is only determined by r2. The 
system enters the regime of unstable quantum modes at 
coupling strength A « ^classic ^ Ac (7). As we shall see, 
this leads to an interesting precursor effect of polaronic 
behavior. 



B. Low-temperature self energy 

In this section we discuss the low temperature self en- 
ergy in detail, in order to clarify the limitations of our 
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theory. For weak to intermediate coupling strengths, 
A < Ac, P{r) is peaked about r = with half-width 
of order T, so an expansion in r is justified if AT/t <C 1, 
i.e. if T is low and the system is not too close to the po- 
laronic instability. It is convenient to rewrite Eq. (^) to 
order ^(A^i^,,;,,! 7^,,;,,!, A^j^^^i.^i T^/t^) by adding and 

subtracting the fc = term in Fn and moving all terms 
to the denominator, so (denoting ^ J drP(r) by ()) 



Qn = 



Cn~r ~ S^'*^'*'''(r) - Aclassical(?') T / {Cn ~ r) 



with 



(36) 



(37) 



At low T and not too close to the polaronic instability 
we may proceed by expanding Q to order r^, averaging, 
and then replacing terms in the denominator. This leads 
to a self energy 



^lowT 



(38) 



The expression for An is somewhat cumbersome, because 
many terms contribute and the result depends on the 
relative values of external and phonon frequency. How- 
ever, in all cases we have examed, the analytic contin- 
uation of An — !■ A{uj + iO^) has an imaginary part of 
positive sign, corresponding to a non-causal contribution 
to S, and magnitude of order unity (7''). For exam- 
ple, if one assumes that all relevant contributions in the 
occuring Matsubara sums stem from small frequencies 
<<2t then (noting c„ = itsgnujn at half filling) 



. / „N . Aclassical , ■ \ \ 

A[UJ = 0) = I ; 1- I Aclassic A 

7classical 



o ■'^classic -'^4 + Adassic A 



(39) 



The non-causal self energy is of order C'(7ciassicai T), 
i.e. formally C'(7ciassicai)' ^^'^ ^^^^ from the deriva- 
tion to arise from an incomplete treatment of the 0(7^) 
contributions to physical quantities. Unfortunately at 
(w, T) <C r2, the leading contribution Siviigdai ^ ^ s{Q/T) 
where s ~ exp[— il/T] as T — > 0, so the non-causal cor- 
rection term dominates. 

This may be viewed in a different way: a purely classi- 
cal approximation (Q —^ 0) would lead to a T-linear scat- 
tering rate with coefficient of order unity. The semiclas- 
sical approximation overcorrects for this behavior, can- 
celling the leading term leaving a correction of order jT 
but of the wrong sign. 

The unphysical low temperature behavior is not due to 
the realization of the semiclassical approach introduced 



in E g. ( pif). Starting from the full expressions Eqs. ( |T5| ) 
and (|l^) we find a similar result 



A, 



si, - A ss A^ 



,A2{2s + l)T^x 



n fc>0 



1 



1 



1 



1 1 

^n+k 



(40) 



which will also lead to a non-causal self energy. We con- 
clude, that higher orders of quantum contributions 7" 
are required to completely suppress the thermal contribu- 
tions at low temperatures T < D, (which wrongly would 
lead to a linear resistivity in T). 

However, the thermal fluctuations are reduced in the 
semiclassical approach compared to the classical (7 = 0) 
case by a factor 7 and the low-temperature self-energy is 

M) and (p9|): 



of order as can be seen from Eqs 



S (w^O) 



. 

classic 

t 



f 



T 



n, 



lassie 



(41) 



where / is a scaling function. Therefore, an order 7ciassic 
impurity scattering rate r, as proposed in the previous 
section, will lead to a well-defined electron self-energy. 
An appropriate choice may be found from Eq. (p8[). In 
the classical limit (ilciassic 0) we find 



(42) 



i.e. a scattering jate r « Adassicai ^^ciassicai/12 for tem- 
peratures T < r^ciassicai- This is the value we choose 
for the additional impurity contribution. Note that in 
the quantum limit (T — > 0) the static approach gives 

T ~ Aclassical ^classical/2. 



C. Observables 

We end this section by introducing the four observables 
we have investigated numerically: phonon distribution 
function, electron spectral function, optical conductivity, 
and resistivity. The normalized phonon distribution func- 
tion pir) = P{r) /Z was defined above. Within DMFT 
the spectral function and optical conductivity can be de- 
rived from the electron self-energy on the real frequency 
axis. The density of states is given by 



p{uj) 



de po(e)p(e,w) 



with the spectral function defined by 



1 



■ Im 



1 



w + ^ - Y,{uj) - eg 



(43) 



(44) 



The optical conductivity is obtained from linear response 
theory applied to the expectation value of the current 
operator 
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^ = E ( ^^^^ ^''^ 

In performing the sums ovKr momentum p we will as- 
sume a hypercubic lattice.E3 The differences to a Bethe 
lattice are not significant for the results presented here. 
The locality of the electron self energy implies that the 
momentum summation can be performed independently 
on each side of an electron-phonon vertex when calculat- 
ing the expectation value of the current. Since the elec- 
tron velocity dtp/ dp is an odd function of momentum all 
vertex corrections vanish and the optical conductivity is 
obtained by a convolution of two full Green functions: 

(T(ri) = (ToX 

where ctq is a constant and /(lu) is the Fermi distribu- 
tion. The resistivity of the system can be deduced from 
the optical conductivity via puc = 1/^(0)- We briefly 
discuss Puc in the regime A < A^^'^*'*'"^ where the quan- 
tum modes are stable. If E = 0) <C 2t as T 0, the 
low temperature DC conductivity is given by 

--W = ^"^w/^ eosh^t/2T) Mb)| ^''^ 

and we expect (for a Fermi liquid ground state) ^ 
const -I- T^. However, due to the low-temperature defect 
of the semiclassical approach discussed above this regime 
cannot be reached. 



V. RESULTS 

In this section we apply the semiclassical method de- 
veloped above to a system of spinless electrons (s = 0) 
interacting with local lattice distortions, which is rele- 
vant to half-metals such as CMR manganites. First we 
consider high to intermediate temperatures. As discussed 
in section III.B the high to low temperature crossover of 
dynamical quantities can be discussed within the 'i^ = 0' 
implementation of the semiclassical approach. This will 
be done first. A detailed comparison with other ap- 
proaches will be given below, when discussing the low 
temperature behavior. 

Fig. ^ shows the resistivity as a function of tempera- 
ture for a wide range of coupling strengths and 7 = 0.1. 
For narrow electron bands (typical for e.g. the A-15 ma- 
terials) this choice corresponds to setting ~ 300K. At 
high temperatures T ::§> Q mu^produce the results of the 
classical (7 = 0) approach.oEj At weak couplings the re- 
sistivity depends linearly on temperature p(T) = AT+B. 
At intermediate couplings the resistivity becomes non- 
linear in T and the ratio A/B decreases. As discussed 
in detail in Ref. this is essential for the phenomenon 
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FIG. 1. Resistivity as a function of temperature T for various 
electron-phonon coupling strengths A at a given phonon frequency 
7 = 0.1. Temperatures corresponding to the renormalized phonon 
frequency |n(r = 0)| (taken at T = 0.2) are indicated by a cross for 
each value of A. Dashed and solid lines correspond to the regime 
of stable (A < A^'''^^) and unstable (A > A^'^^^) quantum modes, 
respectively. The low T behavior in the regime of unstable quantum 
modes is beyond the scope of this paper and the dotted line for 
A = 2.89 serves as guide to the eye. 

of 'resistivity saturation' in models which couple local 
lattice vibrations to the level positions: the scattering 
rate of the electrons with lattice fluctuations increases 
with T, but less rapidly than predicted from second or- 
der perturbation theory. Notice that when discussing re- 
sistivity saturation we focus on temperatures above the 
(renormalized) Debye temperature indicated by a cross 
in Fig. ^ As seen in Fig. |l| the system remains metal- 
lic for coupling strengths well above A^'*^'^*'"^ — 2.36 be- 
cause quantum lattice fluctuations allow the electrons to 
tunnel between neighboring sites. The crossover to the 
high temperature resistivity (which is large) is more pro- 
nounced than for a conventional metal. We thus refer 
to this state as a bad metal. Our approach is capable of 
describing this interesting state where quantum (Migdal) 
and thermal (classical) fluctuations compete. The strong 
renormalization of the phonon frequencies leads to a drop 
in resitivity at a temperature which depends strongly on 
A and which is well below the unrenormalized Debye tem- 
perature 0.1. Only for stronger couphngs A > 3.03 
does the resistivity show insulating behavior over the en- 
tire accessible temperature range. 

The new feature brought by quantum fluctuations to 
these curves is a coexistence of very weak temperature 
dependence at high T with metallic behavior at low T. 
This is clearly shown in Fig. ||, which displays the ef- 
fect of changing quantum fluctuations on the A — 2.89 
curve of Fig. 0. One sees that the nearly classical model 
(7 = 0.01) is an insulator at low T, while the other two 
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FIG. 2. Resistivity as a function of temperature T for interme- 
diate electron-phonon coupling strength A = 2.89 > A^'^*^*^ = 2.36 
and various phonon frequencies. Quantum fluctuations destroy the 
polaronic insulator and restore a metallic state. 



are metals, while the higher-T behavior is hardly affected. 
The 7 = 0.1 data bear a striking resemblance to the resisa 
tivity of well known saturation materials such as NbsSb.Q 
We also note that the data shown in Fig. ^ imply an 
isotope- driven metal insulator transition, although to ob- 
tain this coupling one must finetune A within 0(7) of the 
critical value. 

In Fig. we indicated the modulus of the renormalized 
phonon frequency by a cross. However, |f2(r = 0)| may 
be considerably temperature dependent. Fig. [s] shows 
the temperature dependence of the renormalized phonon 
frequency \n{r = 0)|. For weak couplings (e.g. A = 1) the 
frequency is slightly suppressed and nearly independent 
of temperature. However, close to the (classic) polaronic 
instability A^^''''''''^ (e.g. A = 1.96) the phonon frequency 
is strongly renormalized towards zero and temperature 
dependent. Nonmonotonic temperature dependence of 
phonon frequencies was found in models with interactions 
between phonons and strongly correlated electrons and 
was used to interpret anomalies in Raman scattering and 
acoustic expepipients on certain superconducting molec- 
ular crystals .Ell Here we show that even pure electron- 
phonon coupling may lead to a nonmonotonic temper- 
ature dependence of the phonon frequency. Note, how- 
ever, that the connection between |0(r — 0)| and measur- 
able quantities remains yet to be established. On enter- 
ing the regime of unstable quantum modes the temper- 
ature dependence of the renormalized phonon frequency 
change qualitatively. In addition, an enhancement of the 
phonon frequency is possible. In the insulating regime 
|il(r = O)l(r) increases monotonically. 

The collapse of the polaronic state due to quantum 
fluctuations causes other changes in electronic proper- 
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T 

FIG. 3. Renormalized phonon frequency \fi{r = 0)| as func- 
tion of temperature T for various coupling strengths A at a given 
phonon frequency 7 = 0.1. Close to the (classic) polaronic in- 
stability a clear temperature dependence is visible. For A = 1.96 
the relative change of |f2(r = 0)| amounts to 12% over the plotted 
temperature regime. Solid and dashed lines indicate the regime of 
stable [real fl{r = 0)] or unstable [imaginary Q(r = 0)] quantum 
modes, respectively. 
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FIG. 4. LOWER PANEL: Distribution of local distortions 
^'static!'") for three different phonon frequencies and T = 0.2 (dot- 
ted lines), T = 0.1 (dashed Unes), and T = 0.04 (solid lines). The 
coupling strength equals A = 2.89. Scales on x- and j/-axes are 
the same in all plots and given in the first one. UPPER PANEL: 
Electron spectral function as function of frequency. The parameter 
values are identical to those in the lower panel. Scales on x- and 
j/-axes are the same in all plots and given in the first one. 
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FIG. 5. Optical conductivity fr(tt') for three different phonon 
frequencies and T = 0.2 (dotted lines), T = 0.1 (dashed lines) and 
T = 0.04 (solid lines). The coupling strength equals A = 2.89. 
Scales on x- and y-axes are the same in all plots and given in the 
first one. 



ties which are summarized in Figs. ^ and ^. The lower 
panel of Fig. |^ shows the distribution of local lattice 
distortions for various frequencies and the uppper one 
the spectral function. Fig. ^ displays the optical con- 
ductivity. At high temperatures we find a considerable 
probability for finite (dynamic) lattice distortions even in 
the metallic state. The density of states at the Fermi en- 
ergy is relatively small and the Drude peak in the optical 
conductivity is very broad. In fact, the optical spectral 
weight can be spread over a frequency range which con- 
tinously increase with increasing temperature. Therefore 
the term "saturation" is a misnomer: There is no intriapc 
maximum value of the high-temperature resistivity.E2IE3. 
Upon lowering the temperature we recover a conventional 
metallic state, where the probability of strong lattice dis- 
tortions decreases, the density of states at the Fermi en- 
ergy increases, and the Drude peak sharpens. 

In a polaronic insulator there are pronounced nega- 
tive (positive) distortions indicating that the electrons 
(holes) remain on a given site for a sufficiently long time 
for the lattice to relax. If the phonon frequency increases 
the probability at r w becomes larger, which points 
to an increased electron mobility. As a consequence, 
three maxima in Pstatic(?') may be observed when cross- 
ing over from high to low temperatures for sufficiently 
large phonon frequency. The upper panel of Fig. ^ shows 
the corresponding changes in the spectral function of the 
electrons. The polaronic insulator exhibits a (pseudo) 
gap. The 'splitting' of the band is caused by dynamical 
local distortions. Upon increasing 7 one sees that at low 
T considerable spectral weight is transferred from high 
to low energies, filling up the gap, whereas the spectral 
function at high T is virtually unchanged. Further insight 
may be obtained by considering the optical conductivity 
(Fig. ^). For 7 = 0.01 we see an insulating gap develop 
from a broad, incoherent high-T conductivity; again in- 
creasing the phonon frequency leads to the development 
of a low T metallic state. 
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Resistivity as a function of temperature T for A = 1.96, 
7 = 0.1, calculated using various implementations of the semiclas- 
sical approach: (i) Static approximation Ilfc IIo (static, dotted 
line), (ii) expansion in w^, (dynamic, dashed line), (iii) neglecting 
F (solid line, left panel) and (iv) adding an impurity scattering 
(solid line, right panel). For details see text and section III. The 
deficiencies of approaches (i) (nonmonotonic resisitivity) and (ii) 
(negative resisitivity, i.e. non-causal self energy) are clearly visible. 



The technical difficulties mentioned in the previous 
section limited the calculation shown in Fig. | and |^ 
to relatively high temperatures. We now turn to the low- 
temperature behavior (T < f2) to reveal the strengths and 
weakness of different low-T implementations of the semi- 
classical approach. We focus on (A, 7) — (1.96,0.1). Fig. 
^ shows the low-temperature resistivity, calculated using 
various implementations of the semiclassical approach. 
Quantum lattice fluctuations enter in two places: the 
distribution function of lattice distortions P{r) and the 
generalized contribution F„ (r) to the electron self energy. 
The two properties influence the resistivity in opposite 
ways. A higher probability of small distortions decreases 
the resistivity whereas a negative F"{oj,r) describes ad- 
ditional electron-lattice scattering. 

The implementations discussed here only differ in the 
choice of F as discussed in section lll.B. The static ap- 
proximation gives F„(r = 0) = (Adassic f^classic/2) c,7^ - 

AciassicTc^^ at very low temperatures and leads to a 
nonmonotonic resistivity, as seen in Fig. |^ (dotted line). 
When we take into account that F has to be reduced for 
energies around the Fermi energy, the finite resistivity 
at zero temperature of the static approximation is elimi- 
nated but the resulting electron self energy now violates 
causality. Fig. |^ shows that the resistivity drops to zero 
at a finite temperature (dashed line). We conclude that 
at very low temperatures dynamic corrections of higher 
orders (see section lll.B) become important. To mimic 
their effect we suggested, first, to set F = 0. This is good 
at high to intermediate temperatures but cannot account 



9 



0.36 



0.34 



0.28 



0.24 













0.2 








1 — ' — t — ' — r 






— ' — t — ' — 
















0.1 


o 
11^ 

3 


- 









u 







LI- 


1" _ Q .] 

T = 0.08 












-0.1 




T = 0.06 

T = 0.04 




3 -2 - 


1 1 


2 ; 








CO 




/' 


\ 


T = 0.02 










/ 


\ 
















^ ^ \ 























-0.4 -0.2 0.2 0.4 



CD 

FIG. 7. MAIN PANEL: Density of states as function of 
frequency near the Fermi energy for (A, 7) — (1.96,0.1) and 
various temperatures. The polaronic feature which appears 
at ~ is responsible for the strong decrease of the resistiv- 
ity. INSET: Imaginary part of the generalized contribution 
of the quantum lattice fluctuations to the electron self energy 
F{uj,r) as function of frequency u} at zero distortion r — 
(parameter as in the main panel and T — 0.02). Energy con- 
servation reduce the scattering in a shell around the Fermi 
energy. 

for a Fermi liquid like resisitivity, pDC ^ const + as 
seen in the left panel of Fig. ^ (solid line). Second, we 
suggested adding an impurity scattering of order 7 to ob- 
tain a well defined electron self energy. This is shown in 
the right panel of Fig. ^ (solid line): at very low temper- 
atures the resistivity levels off to a finite value due to the 
added impurity scattering. Both ('F = 0' and impurity) 
approaches demonstrate the development of a resistivity 
p ~ at A = A^'''"'^''^'*' and T ^ 0. 

The inclusion of impurity scattering is an ad hoc rem- 
edy. However, it can be used to demonstrate the for- 
mation of polaronic bands. In the calculation leading to 
the spectral function in Fig. ^ we neglected the gener- 
alized contributions of the quantum lattice fluctuations 
to the electron self energy completely, i.e. F = 0. In a 
bad metal, as seen in Fig. ^, this leads to an increase 
of spectral weight in a broad energy range around the 
Fermi energy when the temperature is lowered. How- 
ever, F should be frequency dependent. Additional scat- 
tering of electrons with quantum lattice fluctuations is 
possible for electrons sufhciently far from the Fermi en- 
ergy as enforced by energy conservation (absorption and 
emission of finite frequency phonons). The frequency 
dependence of the scattering rate is illustrated in the 
inset of Fig. ^ which displays the imaginary part of 
F{oj,r = 0) as given in Eq. (|27|). Even the sign changes 
near lu — which compensates partially scattering from 
thermal fluctuations. This physics is captured by the im- 
purity approach and leads to a clearly visible polaronic 



band in the spectral function as seen in the main panel 
of Fig. 0. For weak to intermediate coupHng strengths, 
A < A^'^^'^'^, the phonon distribution function Pstatic('') 
is concentrated around r « 0. Therefore, at T = the 
electron self energy is given by the Migdal expression, 
see Eq. (^8|) , and a similar polaronic peak in the spectral 
function like here at finite temperatures can be observed 
at zero temperature.E3 However, the three peak struc- 
ture of P(r) as seen in Fig. ^ developing close to Ac (7) 
is clearly beyond the Migdal approach. 



VI. CONCLUSION 

In this paper we have developed a semiclassical ap- 
proach based on the dynamical mean field theory to treat 
the interactions of electrons with local dynamic lattice 
distortions. The method is not restricted to small or 
strong electron-phonon couplings. It can be applied to 
three spatial dimensions and may be extended to in- 
clude interactions other than electron-phonon. The ef- 
fective action is organized in terms of local static n-point 
correlation functions r„ which have to be appropriately 
modified in the presence of other interactions. However, 
the frequency dependence of r„ cannot be neglected all- 
together. At low temperatures it leads to additional 
electron-lattice scattering which is important when cal- 
culating dynamical quantities. We discussed the effect 
of these dynamic corrections extensively. Nevertheless, 
errors of order 7^ inherent in the semiclassical approach 
allow only for a qualitative discussion of the low temper- 
ature behavior by introducing an impurity scattering of 
order 7 by hand. A quantitative analysis is beyond the 
scope of this paper and left for future work. However, at 
a critical couphng A — X'^^^^^'^ — 37r/[(2s-|- 1)4] a classical 
picture should hold and the resistivity becomes p ^ T^/^ 
for very low temperatures. 

We have applied our method to study isotope effects 
on electronic properties in the crossover regime between 
Fermi liquid and polaronic behavior. We found large ef- 
fects including an isotope-driven insulator-to-metal tran- 
sition. The phenomena described in Sec. V may be 
relevant for various half-metals such as CMR mangan- 
ites and A15 compounds even though the simplicity of 
the model studied here does not allow for detailed com- 
parisons with experiments. Nevertheless, an isoi|©pe- 
driven insulator-to-metal transition was observeda in 
(Lao. 25Pi'o. 75)0. TCao.sMnOs. Within the picture pro- 
posed above this can be understood as a quantum tun- 
neling effect. In addition to electron-phonon interac- 
tions, a realistic description of conduction electrons in 
the CMR manganites must account for double exchange 
with the core spins of the Mn ions. This would am- 
plify the insulator-to-metal transition observed above: if 
the electrons become mobile the core spins tend to order 
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ferromagnetically in order to further lower the electron 
kinetic energy. 

We have also addressed the question of resistivity sat- 
uration in metals with high resistivity. Quantum lat- 
tice fluctuations shift the polaronic instability to larger 



couplings S\ 
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Therefore, systems with strong 



quantum lattice fluctuations remain metallic up to larger 
electron-phonon couplings than systems with purely clas- 
sical (thermal) fluctuations which implies a stronger vio- 
lation of the validity of second order perturbation theory. 
We may (misleadingly) speak of resistivity saturation al- 
though there is no upper bound for the high-temperature 
resistivity in contrast to models which couple jEtomic vi- 
brations to electron hopping matrix elements. t3 The re- 
sistivity at high temperatures assumes very large val- 
ues and a pronounced change in slope of the tempera- 
ture dependent resitivity is obtained when crossing over 
from high to low temperatures. The smoothness of the 
crossover can be changed by isotope replacements with- 
out changing the strength of the electron-phonon cou- 
pling. 
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